A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
10(07.03.01, 07.03.02, 13.09.5) 



ASTRONOMY 

AND 
ASTROPHYSICS 



An ultra-deep ISOCAM observation through a cluster-lens * 

B. Altieri 1 , L. Metcalfe 1 , J.P. Kneib 2 , B. McBreen 3 , H. Aussel 4 , A. Biviano 5 , M. Delaney 3 , D. Elbaz 4 , 
K. Leech 1 , L. Lemonon 4 , K. Okumura 6 , R. Pello 2 , and B. Schulz 1 

1 ISO Data Centre, Astrophysics Division, Space Science Department of ESA, Villafranca del Castillo, P.O. Box 50727, E-28080 
Madrid, Spain. 

2 Observatoire Midi-Pyrenees, 14 Av. E. Belin, 31400 Toulouse, France. 

3 Physics Department, University College Dublin, Stillorgan Rd., Dublin 4, Ireland. 

4 DSM/DAPNIA/SAp, CEA-Saclay, 91191 Gif-surYvette Cedex, France. 

5 Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, 1-34131 Trieste, Italy. 

6 Institut d'Astrophysique Spatiale, Bat. 121, Universite Paris-Sud, F-91405 Orsay, France. 



, received October 1998; accepted asap 

Abstract. We present results of ultra-deep ISOCAM ob- 
servations through a cluster-lens at 7/xm and 15/im with 
the Infrared Space Observatory (ISO) satellite. These ob- 
servations reveal a large number of luminous Mid-Infra- 

■ red (MIR) sources. Cross-identification in the optical and 
Near-Infrared (NIR) wavebands shows that about half 
of the 7 fxm sources are cluster galaxies. The other 7/xm 
and almost all 15 /im sources are identified as lensed dis- 
tant galaxies. Thanks to the gravitational amplification 
they constitute the faintest MIR detected sources, allow- 

. ing us to extend the number counts in both the 7 and 
15 /mi bands. In particular, we find that the 15 /im counts 
have a steep slope ais^m = —1.5 ± 0.3 and are large, with 
5pm(> 30/iJy) — 13 i 5 arcmin 2 . These numbers rule 

■ out non-evolutionary models and favour very strong evo- 
lution. Down to our counts limit, we found that the re- 
solved 7 /im and 15 /xm background radiation intensity is 
respectively (2±0.5)xl(T 9 and (5±l)xl(T 9 Wm^sr" 1 . 

Key words: Galaxies: abundances - Galaxies: clusters: 
general - Galaxies: evolution - Gravitational lensing - In- 
frared: galaxies 



1. Introduction 

Great progress in the understanding of physical properties 
of galaxies has been achieved with Mid-Infrared (MIR) 
and Far Infrared (FIR) observations using the ISO satellite 
(Kessler et al., 1996) and its ISOCAM camera (Cesarsky 
et al., 1996). 

Deep optical surveys, such as the Hubble Deep Field 
(Williams et al., 1996) revealed a new population of dis- 

Send offprint requests to: B. Altieri: baltieri@iso.vilspa.esa.es 
* Based on observations with ISO, an ESA project with in- 
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Kingdom) with the participation of ISAS and NASA 



tant sources at high-redshift (z > 2.5) either using 'drop- 
out' techniques (Steidel et al., 1996) or 'photometric red- 
shift' (e.g. Lanzetta et al., 1996, Connolly et al., 1997). 
Yet optical surveys may miss a whole class of high-redshift 
dust-enshrouded galaxies (Blain & Longair 1993). In par- 
ticular the progenitors of today's E/SO galaxies at z >2 
may emit an important fraction of their total light in the 
NIR, MIR, FIR and Sub-millimeter regions (e.g. Frances- 
chini et al., 1994). Indeed, galaxy formation models (e.g. 
Arimoto & Yoshi 1987, Guiderdoni et al., 1997) predict 
that galaxies in their forming phase are heavily obscured 
by abundant gas and dust inside the system. These is- 
sues have motivated very deep ISOCAM observations of 
blank fields (Rowan- Robinson et al., 1997; Aussel et al., 
1997; Aussel et al., 1998; Desert et ah, 1998; Taniguchi et 
al., 1997), as well as Sub-millimeter observations (Smail, 
Ivison and Blain 1997; Blain 1997; Barger et al., 1998.; 
Hughes et al., 1998). 

The ISOCAM-HDF observations (Oliver et al., 1997; 
Aussel et al., 1998) showed that a non-evolving model at 
15/im is ruled out at 3er. Analysing galaxy counts over 
a wide range of sensitivity, Elbaz et al., (1998) noted a 
change of slope in the counts observed at the mJy-level. 
This variation cannot be explained by a simple evolution- 
ary model, but requires stronger star-formation activity 
for galaxies with 0.4 < z < 1. 

Very deep ISOCAM observations have been conducted 
through gravitationally lensing clusters as part of a large 
programme executed in the ISO Science Operations Team 
guaranteed time (the ISO ARCS programme). The clus- 
ters used were A370 (Metcalfe et ah, 1997, 1999), A2218 
(Altieri et al., 1998a; Metcalfe et al., 1999), MS2137, 
C12244 and A2390. Further to a central goal of the ARCS 
programme, fainter detection thresholds were reached in 
a given observation time, thanks to gravitationnal ampli- 
fication. We report here key results obtained for A2390. 

By pushing ISOCAM to its ultimate limits with the 
help of gravitational lensing, we performed ultra-deep ob- 
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Table 1. Log of CAM Observation of the cluster-lens 
A2390. 9 ESA GT data, 4 Lemonon et al., (1998) data. 



servations through the core of probably currently the best 
studied lensing cluster: Abell 2390 (z = 0.23), to obtain a 
magnified view of the background sky. This both increases 
the sensitivity of our MIR maps and reduces the effects 
of source confusion. These new observations supersede in 
area and depth the previous ISOCAM observations of this 
cluster (Lemonon et al., 1998). 

In Section 2 the observations and data reduction tech- 
niques are discussed. Photometry, source identification 
and counts are presented in Section 3. Results are dis- 
cussed in Section 4. Througout this paper, we used fi = 1, 
A = and H = 50h 50 km/s/Mpc. 

2. Observations &: Data Reduction 

2.1. Observations 

The ESA guaranteed time observations reported here were 
allocated one full ISO revolution (science window of 16 
hours), between December 26 and 29, 1997. We used the 
best observational strategy for detection of faint sources, 
based upon in-flight experience. In particular, the field 
was observed in 4 consecutive revolutions in 4 blocks of 4 
hours. The allocated time was equally divided into the 2 
broad-band ISOCAM filters LW2 (5-8. 5^m), centered at 
6.75^m, and LW3 (12-18^m), centered at 15 yum. 

The observational strategy, that reaches the ultimate 
sensitivity (Altieri et al., 1998b), is performed by raster- 
ing the 32x32-ISOCAM detector array in microscanning 
mode. The pixel-field-of-view (PFOV) of 3" per pixel was 
chosen to obtain high spatial resolution at the expense 
of signal, but obtaining a better sampling of the PSF, 
crucial for source cross-identification. The raster step size 
was 7" (2.33 pixels), the minimum value to step out of 
the PSF FWHM in two consecutive pointings. The size of 
rasters was 10 x 10, so that, in the central part of each sin- 
gle raster, 100 different detector pixels sampled each sky 
pixel. Such redundancy was a key factor in ISOCAM deep 
observations. The final maps cover a field size of 2.6' x 2.6'. 

2.2. Data Reduction & Calibration 

The Abell 2390 raw data of both the ARCS and the 
Lemonon et al., (1998) programme were processed to- 
gether following the steps: 

i) Dark subtraction using an IDL based package (Ott et 
al., 1997) that uses a time-dependent dark correction (Bi- 
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Fig. 1. LW2 image of the A2390 cluster. The beam size 
(defined as 80% of encircled energy) is indicated in the 
lower-righthand corner. Objects detected with the SEx- 
tractor software are indicated by ellipses. Contours show 
the iso-exposure-time of 5ks, lOks, 15ks, 20ks and 25ks. 



viano et al., 1998). 

ii) Deglitching, flat-fielding and sky subtraction of each 
raster, was performed using the PRETI Multi-resolution 
Median Transform techniques adapted for ISOCAM data 
analysis (Starck et al., 1997), a multi-scale analysis that 
decomposes the different parts of the signals. 
Hi) Field distortion correction (Aussel et al., 1998) for sky 
projection, with a final pixel size of 1" . 
iv) The five (4+1) rasters were then stacked together us- 
ing the drizzling technique (Fruchter & Hook 1998). The 
rasters were aligned using the centroids of the 4-5 bright- 
est sources in each waveband. 

The central parts of the stacked 7 and 15 /im images 
are shown respectively in Figs. 1 & 2. The noise in the 
central part of the image is 0.1 /iJy arcsec -2 in LW2 and 
0.2 ^Jy arcsec~ 2 in LW3. 

The calibration into /iJy was done using the refined 
in-flight values from Blommaert et al. (1998). We took 
into account the transient behaviour of the ISOCAM pixel 
signals, which do not stabilize at faint fluxes, but respond 
with transient drifts that depend both on the background 
level and the source intensity. We calibrated this response 
(in the 13 readouts per raster position) at 60% in LW2 
and 80% in LW3, at various levels by running the IAS 
transient correction model (Abergel et al., 1996). 
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Fig. 2. As Figure 1 but for the LW3 filter 

3. Analysis & results 

3.1. Detection & Photometry 

Source catalogs from our field in LW2 and LW3 were con- 
structed using the SExtractor package (Bertin & Arnouts, 
1996). The detection algorithm searches for 15 contigu- 
ous pixels, each having a surface brightness exceeding a 
threshold (chosen as ~1.5 a of the sky noise, Table 2), 
after subtracting a smooth background signal and con- 
volving the image with a gaussian filter of the same width 
as the 7 and 15 \xm image PSF. We performed 7" aper- 
ture photometry. Aussel et al. (1998) showed that aper- 
ture photometry is very linear, above 100 /iJy even though 
at the faintest fluxes the remnants of cosmic rays pollute 
the sources with a positive bias. The photometry was cor- 
rected for loss of flux in the PSF wings, by using measure- 
ments of calibration stars under similar conditions (mi- 
croscanning & rastering): 60% of the flux is found within 
a 7" diameter at 7 /im and 48% at 15 (im . 

To assess the contribution of noise to our catalogs we 
ran the detection algorithm on the negative fluctuations 
in the maps, concluding that we have no false detections 
above 60 ^Jy in the central 2'x2' area of the maps in both 
filter bands. 

Finally, to determine the completeness of the sample 
we added a template faint source (scaled-down version of 
a calibration star) to the maps repeatedly at different po- 
sitions in the map and estimated the efficiency of detect- 
ing this source as a function of its flux density. This pro- 
vided a reliable estimate of the visibility of a faint compact 
source in the maps. The estimated 80% and 50% com- 
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Table 2. MIR source counts in the image plane, N: num- 
ber of detected sources, N neg : number of detected sources 
on the inverted image 



pleteness limits of the catalogs derived from these simula- 
tions are listed in Table 2. These figures refer to apparent 
source brightness and are similar to the deepest obser- 
vations published to date at 7 (im on the Lockman Hole 
(Taniguchi et al., 1997, faintest detections around 30/xJy) 
and at 15 /zm on the HDF (faintest sources around 50 
/uJy). However, thanks to the gravitational magnification 
of a factor ^2 to 10, the sources are intrisically the faintest 
MIR sources detected to date. 

3.2. Cluster Contamination 

Thanks to our high-resolution images we have been able to 
unambiguously identify almost all the MIR sources with 
counterparts in deep NIR and optical (HST/WFPC2 and 
ground-based) images. The relative astrometric accuracy 
is found to be better than 1" in both filters. In only a few 
cases we suspect that two sources are blended. There is 
one obvious case in the 7 /im map, where the straight arc 
(Pello et al., 1991) is blended with the nearby elliptical 
galaxy. 

30 sources are detected at 7 ^m . Half of them arc 
easily identified as cluster member galaxies (Pello et al., 
1991, Leborgne et al., 1992, Abraham et al., 1996). The 
5-8.5 /jm emission of the cluster galaxies corresponds to 
4.5-6.9 fiuv restframe emission. For E/S0 galaxies it cor- 
responds mostly to the Rayleigh- Jeans tail of their old 
stellar population as in the Virgo cluster (Boselli et al., 
1998). Two stars are identified, and at least ten sources 
are lensed distant galaxies. These lensed sources are all 
detected at 15 fim . 

At 15 fim , 34 sources are detected in the central 
2.25' x 2.25' field. Only three sources are identified as clus- 
ter members: the cD galaxy (Lemonon et al., 1998, Edge et 
al., 1998) and two star-forming galaxies. Based upon spec- 
troscopic or photometric redshifts, all the other sources are 
identified as faint lensed galaxies. Almost all sources for 
which we have spectroscopic redshifts are background ob- 
jects. Although, we can not rule out some of the targets be- 
ing in the cluster, the probability is very small. The detec- 
tion of almost exclusively background sources in the clus- 
ter images demonstrates that at 15 fxm the cluster-core be- 
comes transparent (as in Sub-mm/UV bands). Therefore 
the key feature is that the cluster-core acts as a natural 
gravitational telescope amplifying the flux of background 
sources, typically by a factor of 2. 
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3.3. Source counts 



By correcting for the lens magnification and surface dilu- 
tion effects, contamination by cluster galaxies, and non- 
uniform sensitivity of our maps, we can derive num- 
ber counts at 15 /zm to compare with blank sky counts 
(e.g. in the Hubble Deep field and Lockman Hole). The 
7 /jm number counts are more difficult to derive due to 
the larger contamination by the cluster and because of 
the small number statistics. 

The number density of sources is high with respect to 
the size of the FWHM (~ 6" diameter at 15 /im ), but the 
PSF is well sampled on the final maps, and its shape can 
be used to separate the sources; only two 15 /jm sources 
lie at the location of pairs of suspected high-z galaxies. 
Our counts are not significantly affected by confusion. 

The occasional blending of the sources has not been 
taken into account, but the surface area occupied by bright 
sources is subtracted for the computation of the surface 
density of the fainter ones (ie. other faint sources could 
be hidden by brighter ones). We used the completeness 
of the detection at 15 jim given in § 3. This correction is 
negligible to the 7 /im counts and was not applied. 

Due to the non-uniform sensitivity of our maps (be- 
cause of observation strategy) and the lensing effect, the 
object density per flux bin was computed using gain- 
dependant surface areas, and only the central 2'x2' area 
was taken into account for the faintest fluxes. 

A detailed lensing model of A2390 has been produced 
by Kneib et al., (1998). The lensing acts in two ways on 
the background population of galaxies: 

i) an amplification of the source brightness, typically by a 
factor of 2, but up to 10 near the caustic lines. 

ii) a surface dilation effect of the area probed, which itself 
depends on the redshift; the space dilation is stronger to- 
wards the centre (core of the cluster) and increases with 
source-plane redshift. 

To estimate these factors we used: the spectroscopic 
redshift for 7 objects (Pello et al., 1991, Bezecourt & 
Soucail 1997), and for the rest we use the best redshift 
estimate obtained with photometric redshift techniques 
(Pello, private communication), and/or lensing inversion 
techniques (Kneib et al., 1998). By analysing the case with 
all background galaxies at a mean redshift z = 1.0, we 
checked the dependence of the results on redshift uncer- 
tainties. 



4. Discussion and Conclusion 

The source counts, corrected for cluster contamination 
and lensing effects, in both the LW2 and LW3 bands, are 
presented in Figs. 3 and 4. We have used only 5cr sources 
(i.e. ~60 /iJy before lensing amplification correction). Ac- 
cording to Hogg & Turner (1998), this is sufficient to avoid 
the positive flux-estimate Eddington bias which occurs for 
faint source counts. 
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Fig. 3. 7 lens-corrected counts, of identified field 
galaxies, (filled stars) compared to the Lockman Hole 
counts (filled squares, Taniguchi et al., 1997). 
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Fig. 4. 15 /jm lens-corrected counts corrected for incom- 
pleteness, full model (filled stars), with all sources at z=l 
(open stars), without lensing correction (filled circles) and 
compared to the HDF counts (Aussel et al., 1998)(filled 
squares) and the non-evolution model (dashed line) from 
Franceschini et al., (1997). Error bars include both Pois- 
son and systematic terms. 
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At 7 /im , we find that the source density is greater 
than that of Taniguchi et al., (1997) by a factor of 2. 
We suggest that their observations are incomplete below 
50 /iJy, whereas our 7 /im map is 80% complete down to 
30 /iJy. 

At 15 /xm the number counts are compatible with the 
results of Aussel et al. (1998) in the HDF. However, with 
the help of gravitational lensing we are able to extend the 
counts in both bands down to 30/iJy. Putting all back- 
ground sources at z = 1 only slightly decreases the faintest 
counts, at the faintest end, because a few suspectedly high- 
z faint sources would be less amplified. This shows the 
small dependency of our derived counts on the models. 
The very small positive magnification bias in the lensing- 
uncorrected counts is consistent with the measured slope 
of the source counts, while the general shift of the points 
towards higher flux follows from the loss of the lensing 
gain. 

We find a total number density Nj(> 30fiJy) — 3.5 
±1 arcmin~ 2 at 7 /im , and ocis = —1.5 ± 0.3, and 
Ni5(>30/iJy) = 13 ±5 arcmin~ 2 at 15 /im . 

The 15 /jm counts show a steadily increasing excess (by 
more than a factor of 10) with respect to the prediction 
of a no-evolution model (dotted line, Franceschini et al., 

1997) . This confirms the steeper count slope below 1 mJy 
found on the Lockman Hole (Elbaz et al., 1998) and is in 
good agreement with the ISO HDF counts (Aussel et al., 

1998) . The counts are a factor 2-3 higher than the bound- 
aries of the counts coming from an early analysis of the 
background fluctuations in the ISOCAM 15 /im map of 
the HDF (Oliver et al., 1997). In particular, we do not 
seem to detect any sign of flattening of the counts at 
the faintest levels, as expected from evolutionary models 
(Franceschini et al., 1997, Oliver et al., 1997). The slope 
stays close to -1.5±0.3 down to 30 /iJy. This source density 
at faint levels favours extreme evolution models, needed 
to fit the counts at brighter fluxes, as shown by Elbaz et 
al. (1998). 

Integrating the number counts over the whole flux 
range and extrapolating the counts we find that respec- 
tively (5 ± 1) x 10" 9 and (2.0 ± 0.5) x 10~ 9 Win" 2 sr" 1 
is emitted at 15 /im and 7 /im . This implies that a larger 
fraction of the UV/optical background is re-radiated in 
the MIR than in the local universe. 

The absolute astrometry of MIR sources is difficult in 
deep surveys because of the lack of any obvious optical- 
NIR counterparts. In the A2390 field, unambiguous cross- 
identification of more than 90% of the sources was possible 
thanks to a large density of sources and a good sampling 
of the PSF. A number of these sources are well correlated 
with faint galaxies in the visible, some of them having very 
red colors in the NIR. 

15 /im ISOCAM deep imaging is a good way to se- 
lect star-forming galaxies and dusty AGNs which are not 
easy to identify in UV/optical surveys. A more detailed 
analysis of the SED of these MIR detected galaxies will 



be necessary to determine the nature of these sources and 
give an estimate of their SFR. Our observations confirm 
that abundant star formation activity occurs in very dusty 
environments at z ~ 1. Great caution must therefore be 
taken to infer global star formation activity based only on 
UV-continuum or optical luminosities of high-z galaxies. 
The global star formation history can be traced fully only 
if the effect of dust is taken into account in a consistent 
way. (e.g. Blain et al., 1998). Further, detailed analysis of 
our selected sample is needed to unveil the nature of these 
MIR galaxies. 
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